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ABSTRACT 

Context. The collapsar model requires rapidly rotating Wolf-Rayet stars as progenitors of long gamma-ray bursts. However, Galactic 
Wolf-Rayet stars rapidly lose angular momentum due to their intense stellar winds. 

Aims. We investigate whether the tidal interaction of a Wolf-Rayet star with a compact object in a binary system can spin up the 
Wolf-Rayet star enough to produce a collapsar. 

Methods. We compute the evolution of close Wolf-Rayet binaries, including tidal angular momentum exchange, differential rotation 
of the Wolf-Rayet star, internal magnetic fields, stellar wind mass loss, and mass transfer. The Wolf-Rayet companion is approximated 
as a point mass. We then employ a population synthesis code to infer the occurrence rates of the various relevant binary evolution 
channels. 

Results. We find that the simple scenario — i.e., the Wolf-Rayet star being tidally spun up and producing a collapsar — does not 
occur at solar metallicity and may only occur with low probability at low metallicity. It is limited by the widening of the binary orbit 
induced by the strong Wolf-Rayet wind or by the radius evolution of the Wolf-Rayet star that most often leads to a binary merger. The 
tidal effects enhance the merger rate of Wolf-Rayet stars with black holes such that it becomes comparable to the occurrence rate of 
long gamma-ray bursts. 
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^ 1. Introduction 

'NT . 

T-H , Although our understanding of long gamma-ray bursts (GRBs) 
' has increased significantly since they were first discovered, it 

C5 I is still not clear what their progenitors exactly are. From stud- 

T^ljj- . ies of host galaxies it emerged that they occur in or near star- 
' forming regions an d that several are associated with energetic 

OO Type Ic supernova dHiorth et al.i .2003). These supernovae are 
thought to stem from to the explosion of a massive Wolf-Rayet 
^ ] (WR) star, although the signatures of a WR type progenitor have 

• I— I . been unam biguously found in the afterglow of only one GRB 

^ ■ (Ivan Marie et al.ll2005l) . 

JH ' The most widel y used model fo r the formation of GRBs is 
. . the collapsar model (IWooslevlll993h . In this model the core of a 
massive, fast rotating star collapses into a black hole. An ac- 
cretion disk is formed around the black hole if the core has 
enough specific angular momentum, i.e. 7 > 3 x lO'^cm^s ' 
jMacFadven & Wooslevl[T99ft) . The remainder of the core is ac- 
creted onto the black hole and a highly relativistic collimated 
outflow is produced which releases a large amount of energy 
(lO*"' ergs). If the star has no hydrogen envelope, the light cross- 
ing time is less or comparable to the duration of the accretion. 
In this case a GRB can be formed along with a Type Ib/c super- 
nova. Thus, rapidly rotating WR stars are required to produce a 

collapsar. 

Earlier work on GRB progenitors (iPetrovic et al.l 120051: 
iHirschi et al.l2006l) has shown that stellar models without a mag- 
netic field can have enough specific angular momentum in their 
cores to produce a GRB within the collapsar model. The single 
star models in which magnetic fields according to Spruit (2002) 
were included had too little specific angular momentum in their 
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cores due to the increased cor e-envelope coupling , which leads 
to a spin-down of the core (iPetrovic et al.1 l2005t iHeger et aP 
2005). The reason why magnetic fields are considered to be im- 
portant are the low observed rotation rates of young neutron stars 
and white dwarfs. In order to rep roduce those low rates, mag- 
netic t orque need to be inclu ded dHeger et al.ll2005t ISuijs et aP 
120081) . IPetrovic et al.l (l2005b also considered binary models in 
which the secondary star was spun up to close to critical rota- 
tion due to mass transfer from the primary. In their model, ini- 
tially the core of the companion was spun up due to magnetic 
core-envelope coupling, but the same mechanism decreased the 
core angular momentum by almost a factor of 100 once the star 
had reached core helium burning. Therefore, either one has to 
consider stars with a low metallicity which have a lower mass 
loss rate, or GRBs at solar metallicity require a more exotic bi- 
nary channel. The first opti on has been considered recently by 
lYoon & Langeii (l2005h and lYoon et al.1 (|2006|) and the second 
option is discussed here. 

We investigate here whether tides in a close binary system 
can spin up a WR star enough so that it can form a collapsar The 
magnetic fields inside the star enforce close to rigid rotation dur- 
ing core helium burning, so that angular momentum added to the 
envelope due to the tidal interaction is expected to be transported 
to the core of the WR star. In order to have a sufficiently strong 
tidal interaction for spin-up to occur, the orbital period needs 
to be smaller than about 24 h. This restricts potential compan- 
ions to compact stars, helium stars, or low mass main-sequence 
stars. While the number of double- WR star systems is too small 
to be significant ( V rancken et al.llT991l) . lvan den Heuvel & YoonI 
(2007) showed that for main sequence companions, the obtained 
spin-up is insufficient to produce a GRB. Therefore, we focus on 
compact companions in the following sections. 
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Currently, only one compact WR binary is known in our 
Galaxy, namely Cyg X-3. The system has a period of 4.8 h 
and a period derivative in the range of 1 - 2 x 10"^ yr ' 
jLommen et a l. 2005). Cyg X- 3 has l ong been known as an X- 
ray binary. Wan Kerkwiik et alj (1 19961) found that the primary is 
a WR-star, probably of the WN spectral type. Th e exact mass of 
the WR star is unknown, but lStark & Saial (12003 ) have placed an 
upper limit on the mass of 7.3 M© using observations taken with 
Chandra. Although the na ture of the companion is still not clear, 
IVan Kerkwiik et alJ (Il996h assume it is a neutron star of 1 .4 M©, 
while Stark & Sai;^ (|2003|) place an upper limit of 3.6 Mq on the 
compact object mass. 

We assume a mass of 10 Mq for the WR star and 1.4 Mq 
for the compact object in our Cyg X-3 case study. The mass loss 
rate is highly uncertain, and estimated by variou s authors to be 
between 0.5 X 10"'' Moyr"' -2.9x10""^ Mq yr"^ dLommen et alJ 
I2005h . 

Systems similar to Cyg X-3, and those considered further 
on in this paper, can only be formed through common-envelope 
evolution. We use Cyg X-3 as a case-study for the tidal spin-up 
process to determine the effect of this process and to see whether 
it is capable of spinning up the WR star enough to meet the col- 
lapsar criterium. The remainder of this paper is organised as fol- 
lows. In Sect. 2 we explain our physical assumptions as well 
as our numerical methods. Our binary evolution models are de- 
cribed in Sect. 3, while Sect. 4 contains our population synthesis 
results. Finally, our conclusions are discussed in Sect. 5. 

2. Methods and physical assumptions 

Our stellar mod els are calculated wi th a hydrodynamic stellar 
evolution code dPetrovic et al.l l2005b . Magnetic fields and the 
transport of angular momentum due to magnetic torque are in- 
cluded ( Spruit 2002), as well as the effects of the centrifugal 
force on the stellar structure, chemical mixing and transport of 
angular momentum due to rotationally induced hydrodynamic 
instabilities, and enhanced mass loss due to close-to-critical ro- 
tation. The stellar wind ma ss loss of the WR star is calculated 
according to H amann et alJ (^.1995, labelled WRO from now on) 
for WR stars with log L/Lq > 4.5: 

log(MwR/M0 yr-i) = ex (-11.95 + l.5\ogL/LQ), (1) 

and lHamann et alj (|1982|) for helium stars with log L/Lq < 4.5: 

log(M„,/M0 yr-i) = e x (-35.8 + 6.8 logL/L^), (2) 

with MwR the mass loss rate of the Wolf-Rayet star and e is a 
multiplication factor used to modify the mass loss rate. 

When a star appoaches critical rotation, the mass loss of the 
star will be enhanced such that over-critical ro tation is prevented. 
To achieve this, we follow the prescription of lLangeii(ll997h : 
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Here V^oi is the rotational velocity of the star, M is the mass loss 
rate of the star, F = kLKAjicGM) is the Eddington factor, and k 
is the opacity coefficient. We apply three different rates of mass 
loss, both based on Eqs. (1) and (2), using e - 0.1, e - 1/3, 



and £ — I. We call these mass loss rates WRl, WR2, and WR3, 
respectively. The WR2 rate roughly agrees with the mass loss 
rate for Galactic WR stars (Nugis & Lamers 2000). We chose 
the first rate to study the effect of a lower mass loss rate on the 
evolution of the system. This essentially captures the situation of 
a lower initial metallicity, sinc e the mass loss rate of a WR star 
is dependant on its metallicity ("V ink & de Koteil2005l) . while all 
other properties depend very little on it. 

For a given WR mass loss rate, one can define an angular 
momentum loss timescale, i.e. the time it takes for the star to 
lose most of its angular momentum to the stellar wind. 
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where Mwr is the mass of the Wolf-Rayet star, Jwr is the total 
spin angular momentum of the star and a ^ 0.1 for efficient 
internal angular-momentum transport, since angular momentum 
loss from the wind is ab out 10 times faster than the mass loss 
(iPacke jl 1 98 ll: iLangeJ 1 998.) . 

There are two different mechanisms to synchronize a star 
with the orbit, the equilibrium tide and the dynamical tide. 
Although both mechanisms create a tidal bulge which causes 
the binary to become synchronized, they do so on different 
timescales. The time it takes for a star to synchronize with the 
orbit if it has a convective envelope is defined as dZahnll 19771) : 



' sync 



-lid 
q — 
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where q = M^c/Mwr is the mass ratio of the binary, with Mwr 
the total mass of the WR star and Mcc the mass of the compact 
companion, d the orbital separation and R the radius of the star 
to be synchronized. We choose to use this ansatz for the equi- 
librium tide for determining Tsync, although it was originally de- 
signed for convective stars, for two reasons. First, in the radiative 
envelopes of rapidly rotating stars, as considered here, rotation is 
supposed to produce a high level of turbulence, implying that tur- 
bulent viscosity dom inates over the radiative viscosity. Second, 
iToledano et al.l (1200 7). applying a kinematic model of tidal inter- 
action to early type stars, found that the tidal time scale in these 
models indeed scales as (d/R)^. 

The angular momentum exchange betwe en the orbit and the 
WR spin is computed in the following way dWellsteinll200ll) . If 
AJ is the amount of spin angular momentum that is added to the 
star in the considered time step Af due to tides then. 



AJ - (7wR - -/sync)(l 



0, 



(7) 



where 7sync is the spin angular momentum the star would have 
when the system is synchronized. 

Angular momentum AJ is added to the outer layers at every 
time step. Angular momentum transport processes (i.e. magnetic 
torque) redistribute the angular momentum over the whole star 
Magnetic torque is strong enough to keep the star close to rigid 
rotation during core helium burning. In this way the whole of the 
WR star is spun up due to the tidal interaction. 

To determine the outcome of the common-envelope evolu- 
tion we calculate the final orbital separation and compare that to 
the CO-core radius and the Roche-lobe radius of the CO-core. 
We do this using the energy equation for common envelope evo- 
lution dWebbinkll 1 98l Ide Koollll990h : 
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Table 1. Table 1. Cyg X-3 models. 
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where Of is the final orbital separation after the CE phase 
and a; the orbital separation at the moment when Roche Lobe 
overflow (RLOF) starts. Mcoie is the mass of the CO-core of the 
WR star and Mgnv the mass contained in the envelope of the WR 
star The efliciency parameter ?/ce is set to 1 in our calculations 
and A is determined by the binding energy of the envelope, which 
is extracted from our models. 

3. Binary evolution models 

We first compute tailored models for Cyg X-3 to determine the 
role of tidal interaction in this system. We assume a compact 
object mass of 1.4 (i.e. take a neutron star), chose a WR- 
star mass of 10 Mq and set the initial orbital period to 4.8 hrs. 
The mass loss rates for our Cyg X-3 models are based on the 
iHamann et"an d 19951) rate, but with e - 1 (model WRa), e - 
0.1 (WRb) and e = 0.01 (WRc, see Tableffli. Their evolution is 
described in the next subsection. 

3.1. Cyg X-3 models: a Case Study 

The three Cyg X-3 models which we have calculated show the 
eflBCts of diff'erent mass loss rates on the tidal interaction and 
evolution of the system. The model sequences end either at car- 
bon depletion in the core or when the primary fills its Roche- 
lobe. In all three cases the tidal interaction transfers angular mo- 
mentum from the orbit to the WR star, but the WR star is only 
spun up if the mass loss rate is not too high, i.e. if Tsync < T"ang. 
This is the case for models WRb and WRc, but for model WRa 
the mass loss rate is too high and thus Tang too small. 

Figure [T] shows the mass loss rates for models WRa, WRb 
and WRc as function of time. The WRa mass loss rate de- 
creases throughout the evolution and only increases sharply at 
9.5 X 10^ yr at core helium exhaustion, because the luminosity 
increases greatly when the WR star starts its He-shell burning. 
The WRb and WRc mass loss rates are initially so low that the 
tides spin up the stars to critical rotation, which occurs here well 
before the WR stars fill their Roche lobes since their luminosi- 
ties are near their Eddington luminosities (cf. Eq. |4|. The high 
luminosity of the WR star has the consequence that the star can 
approach critical rotation for wwr < Worb- Once at critical ro- 
tation, the mass loss rate is determine d by the mech anical con- 
straint to avoid over-critical rotation ( lLangeij |l998'. cf.). While 
this is a nume rically (and perhaps physically) unstable situation 
(lLangeiill998 !, cf. Fig. 3 in), the corresponding large oscillations 
of the mass loss rate have no consequences for the evolution of 
the system, as long as the long-term average of the mass loss 
rate is still well defined, as is the case here (cf. Sect. l3.1.2l i. As 
the orbital period decreases, the time-averaged mass loss rates of 
systems WRb and WRc increase in the course of evolution. 

The time evolution of the orbital separation of these three 
sequences is shown in Figure|2] The orbital separation of models 
WRb and WRc decreases, while it increases for model WRa. The 
evolution of model WRa is dominated by its mass loss. The two 
kinks which occur at 3 x 10^ and 9.6 x 10^ yr can be traced back 
to the changes in the mass loss rate at those times (see Figure[T]). 
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Fig. 1. mass loss rate as a function of time for the three Cyg X- 
3 models: WRa (blue line, solid) with e - 1, WRb (green line, 
dotted) with e = 0.1 and WRc (red line,striped) with e = 0.01 
(cf. Table 1). 
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Fig. 2. Orbital separation as a function of time for our 3 diiTer- 
ent Cyg X-3 models. Plotted are model WRa (solid line), model 
WRb (dotted line) and the model WRc (dashed line). 



In models WRb and WRc the initial mass loss rate is low, so that 
the tidal interaction spins the star up and angular momentum is 
thus transferred from the orbit to the WR star. That means that 
the orbit shrinks and this process continues until the WR star 
starts to fill its Roche-lobe. 

3.1.1. Sequence WRa 

Model WRa reaches carbon depletion in the core, although the 
mass of the WR star is only 2.2 Mq at that time. The specific 
angular momentum of model WRa does not change significantly 
during the core helium burning phase, implying that the wind in- 
duced loss of spin angular momentum is compensated by tidally 
induced angular momentum gain. At the end of core carbon 
burning, the specific angular momentum drops to very low val- 
ues in the core. This is due to the overwhelming effect of the 
increased mass loss, which brings the WR star out of tidal lock- 
ing. 
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Fig. 3. Kippenhahn diagram of the specific angular momentum 
for the WR star in model WRa. The hatched and crossed areas 
indicate convection and semi-convection respectively. The red 
color indicates the amount of specific angular momentum, where 
a darker red indicates larger specific angular momentum. 

The internal specific angular momentum distribution during 
the evolution of the system can be seen in Figure [3] The angu- 
lar momentum in the core stays roughly constant throughout the 
evolution, except for an initial increase during the first 40 000 yr 
due to tidal spin-up, and a large decrease after helium shell burn- 
ing has started at 9.5 x 10^ yr When 50% helium is left in the 
core, the star has already lost more than half of its mass due to 
its strong stellar wind. The WR star will not form a black hole 
because the CO core mass is too low (below 2 M©), but it will 
most likely form a neutron star 

3.1 .2. Sequences WRb and WRc 

Model WRb has an entirely different evolution. The star is spun 
up due to tidal interaction and as a consequence its specific an- 
gular momentum increases. At the end of the calculation, the 
specific angular momentum at 3 Mq is lO'^cm^^s"'. Figure |4] 
shows the internal specific angular momentum distribution dur- 
ing the evolution. Because the orbital separation decreases with 
time, the star remains tidally locked and is spun up throughout 
the whole evolution. Figure |4] also shows that despite the oscil- 
lations of the mass loss rate as function of time (Fig. 1), the time 
average mass loss rate is well defined, and in fact dictated by the 
need to avoid faster than critical rotation at the stellar surface. 
Even though the star loses about 3.5 the star is still spinning 
rapidly at the moment RLOF occurs. The outcome of the RLOF 
phase will most likely be a merger, since the mass transfer pro- 
cess is unstable because of the large mass ratio. 

The evolution of model WRc is similar to that of model 
WRb, except that the orbital shrinking is less severe in this 
model. Although this may seem strange at first because the mass 
loss rate is lower and thus the tidal interaction should be more 
effective in spinning up the star, the stronger mass loss in model 
WRb actually helps to decrease the orbital separation. This is 
because the system is tidally locked in both models such that 
any angular-momentum loss from the star comes from the orbit 
instead. This effect can be compared to magnetic braking in low- 
mass stars with a convective envelope. Since the mass loss rate is 
higher in model WRb, the angular momentum loss is also larger 
and the orbit shrinks faster. 
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5.0x10^ 1.0x10= 1.5xi05 2.0x10^ 2.5x10^ 3.0x10= 
t/yr 

Fig. 4. As Figure[3]except this is for model WRb. The transport 
of angular momentum from the envelope to the core can clearly 
be seen after 80000 yr. 
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Fig. 5. As Figure [3] except this is for model WRc. 



Figure |5] shows the internal specific angular momentum as 
function of time. Again, the increase of specific angular momen- 
tum can be seen by the increase of the dark red area in the star. 
Interestingly, the final mass of the WR star is smaller in model 
WRc than in model WRb. This is due to several effects. First, al- 
though the initial mass loss rate is different for models WRb and 
WRc, the effect of the spin-up is that the WR star approaches 
critical rotation. The mass loss rate increases according to Eq.[3] 
The initial value of the mass loss rate does not matter once the 
star approaches critical velocity, which is why models WRb and 
WRc have an almost equal average mass loss rate. Secondly, the 
orbital decay of model WRc is less severe and so the point at 
which the WR star fills its Roche-lobe is reached at a later time, 
so the mass of the WR star is smaller when RLOF starts. The 
likely outcome of this RLOF is again a merger due to unstable 
mass transfer, as is the case for model WRb. 

3.1.3. Summary 

In the three Cyg X-3 models discussed above, either the spin- 
up works and the orbit shrinks (models WRb and WRc) or the 
mass loss rate is too high and the orbit widens (model WRa). 
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The WR star is initially spun up by the tidal interaction in model 
WRa, but after 50000 yr the star starts to lose angular momentum 
due to the strong stellar wind. Models WRb and WRc also spin- 
up initially, but retain their angular momentum throughout their 
evolution. The angular momentum of the orbit always decreases, 
due to the mass loss and the tidal interaction. All three systems 
synchronize their orbits very quickly and remain synchronized 
throughout. 

From these three models it is clear that there is a threshold 
value for the mass loss rate which determines the outcome of the 
system. Either the mass loss rate is too high so the WR star spins 
down and the orbit widens, or the rate is low enough so that the 
WR star is spun up. The mass loss rate is determined not by the 
properties of the WR star, but rather by the need of sufficient 
angular-momentum loss to not exceed the O-limit. Considering 
the limiting case of the mass loss rate approaching zero, it is 
clear that the tidal interaction is strong enough to spin up the star 
and thus the orbit would shrink. The star would then eventually 
approach critical rotation, and the mass loss is enhanced to a 
value which is high enough to keep the star from reaching critical 
rotation. So the spin-up works as shown in the WRb and WRc 
systems, but both systems merge in the end. 

3.2. Parameter study 

It is clear that a study of only Cyg X-3 is insufficient to inves- 
tigate the whole parameter space. We therefore set up a grid of 
binary models with a mass of the WR star between 6 and 18 Mq. 
Our models all have solar metallicity, and we take a zero-age 
helium main sequence star (ZAHeMS) as our WR star We use 
different masses for the compact object: 1.4 M©, 3.0 Mq and 5.0 
M0, i.e. assume either a neutron-star or a black-hole compan- 
ion. While larger compact object masses are possible, systems 
with such are expected to contribute at most 10 - 20 % to the 
total number of helium stars plus compact object systems. The 
initial orbital separation is chosen such that the synchronization 
timescale rsync is about equal to the stellar wind induced angular 
momentum loss timescale Tang. The other chosen initial separa- 
tions are twice and thrice this equilibrium value. For the 5 Mq 
companion we chose an extra set of models, with an initial or- 
bital separation of 0.5 times the equilibrium value. 

As shown for the Cyg X-3 models above, if the spin-up 
works then the orbital shrinking will most likely lead to RLOF 
To investigate whether we can have tidal interaction to spin-up 
the WR star without ensuring RLOF we expanded our calcula- 
tions to several model grids. We use the two different mass loss 
rates labelled as WRl and WR2 (see section|2]i. Each of the plots 
in Fig |6] shows the mass of the WR star and the initial period 
for each system. We computed the evolution of each system and 
gave each different outcome a separate symbol in the plot. We 
found 5 different outcomes for the binary system. 

The first and easiest to understand is the case where the 
initial orbital period is large. The tidal interaction is weak in 
these systems and mass loss dominates the evolution. The binary 
widens and the WR star spins down. Depending on the mass of 
the WR star it will form either a black hole or a neutron star, but 
in neither case is it spinning fast enough to be considered a GRB 
progenitor. 

The second and third type of evolution occur when the tidal 
interaction is not strong enough to spin-up the WR star, so mass 
loss still widens the orbit. However, the system remains com- 
pact enough for the WR star to fill its Roche-lobe during He- 
shell burning. Unstable mass transfer starts at that moment and 
the result is a common-envelope phase (the second one in the 



evolution of the binary). We have found two possible outcomes 
for the common-envelope phase, either a merger, or the binary 
survives the common-envelope phase and ends up as a CO core 
with a compact object in a very small orbit. 

The fourth type is the case where the tidal interaction is 
strong enough to spin up the WR star The orbits shrinks and the 
WR star fills its Roche-lobe during core He-burning. If the mass 
ratio is too far from unity, the mass transfer is unstable and the 
result is a merger of both stars. When the compact companion 
is a neutron star the resu lt will be a Thorne-Zytkow like object 
(iThorne & Zvtkowlll975b . In the case of a black hole compan- 
ion, the WR star will be accreted onto the black hole and a GRB 
may occur (we will discuss this option in the next section). 

In the last type of evolution, the tidal interaction is initially 
strong enough to spin up the WR star to close to critical rotation 
and the orbit shrinks. Due to the WR star approaching critical 
rotation, the mass loss rate increases, which reduces the radius 
of the WR star, resulting in a weaker tidal interaction. Eventually 
the system is no longer tidally locked. The orbit widens again 
and the star continues to lose a large amount of mass. The end 
result is the same as for the systems with a large initial period. 

We do not make a full grid calculation for the WR2 mass 
loss rate models, as the main effect of a higher mass loss rate is 
a shift of the borderlines between the various types of evolution 
in the WRl grids to lower initial periods. The WR2 models are 
similar to the WRl models, but the mass loss is higher and thus 
the initial equilibrium period is smaller. This is because Tang is 
smaller, which means that Tsync also has to be smaller in order to 
have an equilibrium situation. The only potential collapsar pro- 
genitors are those systems which survive the common-envelope 
phase as a CO star with a compact companion in a close orbit, 
if the CO star is spinning fast enough and is massive enough to 
form a black hole. The CO core mass of a 10 WR star is 
6.0 M0, which may be close to the borderline value for form- 
ing a black hole. We have one possible system in our grid for 
which this could be the case, namely the 10 WR star with 
a 3.0 M0 companion and an inital period of 23.15 h (see Table 
lA.U . Out of the 70 systems we calculated with the WRl mass 
loss rate, 19 systems may allow GRB production. The most im- 
portant conclusion drawn from these models is that there is no 
model in which only tidal spin-up works, whilst avoiding RLOF 
or a merger. So our initial idea of having only tidal spin-up is 
unlikely to happen at these mass loss rates and thus at solar 
metallicity. Tidal interaction is very important for triggering the 
RLOF or a merger, thus indirectly contributing to the formation 
of a possible GRB. If tidal interaction is weaker than we have 
assumed here, then the equilibrium period will shift to a smaller 
initial orbital period, the same effect as a higher mass loss rate. 
If the tidal interaction is stronger than we have assumed, then 
the initial orbital period will shift upwards, the same effect as 
a lower mass loss rate. We did not change the strength of the 
tidal interaction, but making the tidal interaction stronger could 
be one way of making tidal spin-up work at solar metallicity, al- 
though the increase would have to be significant because of the 
high mass loss rates at solar metallicity. 

We have investigated the possibility of tidal spin-up in a 
close WR-compact object binary leading to a collapsar. Our re- 
sults show that it is very hard to have tidal spin-up while also 
avoiding a RLOF or merger event. At the same time however the 
tidal interaction is the cause of the pre He-shell burning merg- 
ers. These mergers may or may not produce a GRB, which is 
unknown at the moment. So indirectly tidal interaction may lead 
to the formation of a GRB in these cases. Even if these sources 
would not produce a GRB, they may appear as a transient source 
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Fig. 6. The outcome of the evolutionary calculations for helium star-compact object binaries using the WRl mass loss rate (left 
column) and the WR2 mass loss rate (right column). The companion mass is shown in the upper left corner of each plot. We take a 
1.4 Mq companion (neutron star), a 3.0 Mq companion (black hole) and a 5.0 M© companion (black hole). Evolution outcomes are: 
no RLOF and no spin-up (squares), spiral-in and RLOF during helium burning (crosses), RLOF during He-shell burning resulting in 
a merger (stars), RLOF during He-shell burning resulting in a CO-core + compact object in a close binary (diamonds), and spin-up 
then increased mass loss and orbital widening (triangles). The solid line indicates where the WR star would immediately fill its 
Roche-lobe, which gives the lower boundary for possible initial periods. 
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Table 2. Table 2. Formation rates for each possible GRB pro- 
genitor type, for A - 0.5. 



scenario 


type 


fate 


birthrate [yr ' ] 


A 


He-shell RLOF 


CO-BH merger 


5.64 X lO^" 


B 


pre He-shell RLOF 


He-BH merger 


3.83 X 10-5 


C 


CO + BH 


collapsar? 


1.39 X 10-'' 



on the sky. So our results may have significant observable impli- 
cations, given the fact that our community is currently undertak- 
ing big efforts to investigate transient sources at various wave- 
lengths. 

At lower metallicity (and lower mass loss rates), the range in 
periods at which tidal interaction is strong enough to spin up the 
WR star is larger This is because the 'equilibrium period' shifts 
to a larger initial period, because r^ng is larger and thus also Tjync 
can be larger. There is a limit to which lowering the mass loss 
rate makes a difference, because if the WR star approaches criti- 
cal rotation due to the spin-up, the mass loss also increases and is 
determined by the angular momentum loss needed to avoid ex- 
ceeding the critical rotation rate. So the inital mass loss rate may 
only be important in the initial evolution of the system. One way 
in which there could be only tidal spin-up would be if the initial 
mass loss rate is low and the initial orbital period is large enough 
so that the WR star does not approach critical rotation until the 
end of its evolution. In that way the WR star would have no time 
to spin down again and it will have enough angular momentum 
in its core to form a collapsar. So at lower mass loss rates, having 
only tidal spin-up is still a fine-tuning process, which means that 
it is unlikely that this scenario will be a major GRB production 
channel. Lowering the mass loss rate, i.e. assuming lower metal- 
licity, will have the effect that tidal interaction leading to a RLOF 
or merger occurs more frequently, since the range of periods in 
which this occurs is wider than at solar metallicity. 

4. Population synthesis 

To estimate the birth rate of these kind of systems, we have 
performed a popu lation synthesis study using the model of 
iHurlev et alj (12002.) . A 100^ grid in Mj, M2 and dinn with 
grid limits 13 < Mi/Mq < 55, 6.5 < M2/M0 < 40 and 
60 < (iinit/R© < 1370, which were chosen after searching a 
wide parameter space for progenitor systems, was searched for 
systems which match our GRB progenitors. W e used the ini- 
tial mass function of i Kroupa. Tout. & Gilmor3 (1 19931) for the 
primary, a flat distribution in q = M1/M2 between and 
1 and a separation distribution flat in logc/init between 3 and 
10'* Rq, a Maxwellian supernova kick velocity with dispersion 
190 km/sec, solar m etallicity (Z - 0.02), mass loss according to 
iHurlev et al.l (|2002|) . circular orbits, common envelope param- 
eters Q- = 1.0 and A - 0.5 or 0.05, compact object (post-SN) 
masses according to Belczynski et al. (120021) and the star forma- 
tion rate given in Hurley et al.. (.2002.) of one binary with Mi > 
0.8 M0 per year. 

The formation rates of each of our potential GRB progeni- 
tors are given in Table |2] We only consider systems that have a 
black hole as a companion. These systems can be grouped into 
three categories, each with a different evolutionary path. Group 
A consists of systems that avoid a RLOF phase during the core 
helium burning phase of the WR-star evolution, but experience 
RLOF and a merger during the consecutive expansion of the en- 
velope when helium shell burning starts. Group B contains the 
systems which have a strong tidal interaction during core helium 



burning so that RLOF begins prior to the helium shell burning 
expansion of the star. These systems all end up as mergers. We 
also have one scenario. Group C, which we consider a possi- 
ble collapsar progenitor scenario: The systems that survives the 
RLOF during helium shell burning and end up as CO star with 
a black hole companion. Here, the CO star is likely spun up by 
the tides and the remaining life time may be too short for mass 
loss or expansion of the WR star to be significant. 

The systems in Group B have a significantly higher forma- 
tion rate than systems in the the other two groups. A compar- 
ison of the formation rate of these systems to the average for- 
mati on rate of GRBs pe r galaxy, which is lO-*" yr ' according 
to Cherepashchuk & Post nov (2001), shows that these systems 
could account for a significant fraction of the GRBs or transient 
sources. 

As with all population synthesis studies our results suffer 
from uncertainties in the input distributions and physics. Of par- 
ticular note for the cuiTent work is the highly uncertain common- 
envelope evolution and associated free parameter combination 
aA. We chose those values to match previous studies, such as 
Hurley et al. (2002), but we vary A to show its considerable ef- 
fect. A higher value for A leads to a less tightly bound common 
envelope and a higher chance of envelope ejection — thus en- 
hancing the chance to form close WR star plus compact object 
binaries. Reducing A to 0.05 reduces the rate by a factor of about 
10. This can be understood since a lower A means a more tightly 
bound envelope, i.e. less systems are able to eject the envelope in 
the CE phase and merge. The compact object mass distribution 
(NS/BH) is also quite uncertain but we have no better prescrip- 
tion than that used here. 

5. Conclusions 

This paper shows that it is not a simple matter to spin up a WR 
star through tidal interaction and thereby produce a long gamma- 
ray burst. While the spin-up process itself may work, either the 
evolution of the binary orbit or of the radius of the WR star pre- 
vent the desired result, in almost all cases. Only at low metallic- 
ity, where the WR star winds may be weak and the orbit can thus 
be more stable, can this scenario work for a limited and rather 
insignificant fraction of the parameter space. 

This negative result does not exclude something interesting 
happening to WR stars with a close companion: most systems 
with orbital periods below ~ 20 h lead to a merger In princi- 
ple, the companion could be a compact object, a helium star, 
or a main sequence star. As the spin-up scenario fails, only the 
latter seems interesting in the context of gamma-ray burst for- 
mation (although see Fryer & Heger 2005), in particular when 
the compact object is a black hole. While the merger of a helium 
star with a black hole or neutron star does not form a collapsar, 
it has been proposed th at such events do produce long gamma- 
ray bursts dFrver et al.| [T999). However, detailed models of such 
mergers which can demonstrate their ability to produce a GRB 
are still missing. From the binary models presented above, we 
would expect two types of these events, i.e. a merger during or 
after core helium buring. Both may have rather long time scales 
compared to the average time scale of long gamma-ray bursts. 
It is also unclear whether these mergers can produce an explo- 
sive event resembling a Type Ib/c supernova. Even if these merg- 
ers do not produce a GRB, they may be observed as a transient 
source in the sky. 

In summary, the main product of close WR binaries with 
compact companions is a helium star-compact object merger — 
not a collapsing and rapidly rotating WR star The occurance 
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rate of these events may be compatible with that of long / soft 
gamma-ray bursts. 
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Table A.l. Table 3. Binary model properties for the WRl mass 
loss rate. 



MwR [Mo] 


Mcc [Mo] 








P4 [h] 


6 


1.4 


7.83 


22.16 


40.70 


- 


8 


1.4 


6.53 


18.47 


33.94 


- 


10 


1.4 


5,71 


16.16 


29.68 


- 


12 


1.4 


5.14 


14.55 


26.73 


- 


14 


1.4 


4.72 


13.36 


24.55 


- 


16 


1.4 


4.40 


12.45 


22.86 


- 


18 


1.4 


4.14 


11.72 


21.52 


- 


6 


3.0 


10.40 


29.41 


53.99 


- 


8 


3.0 


8.84 


25.00 


45.93 


- 


10 


3.0 


7.83 


23.15 


40.69 


- 


12 


3.0 


7.12 


20.13 


36.98 


- 


14 


3.0 


6.58 


18.61 


34.20 


- 


16 


3.0 


6.16 


17.44 


32.03 




18 


3.0 


5.82 


16.48 


30.28 




6 


5.0 


12.14 


34.34 


63.08 


4.29 


8 


5.0 


10.49 


29.69 


54.55 


3.71 


10 


5.0 


9.41 


26.62 


48.91 


3.33 


12 


5.0 


8.63 


24.41 


44.85 


3.05 


14 


5.0 


8.04 


22.74 


41.77 


2.84 


16 


5.0 


7.57 


21.41 


39.34 


2.68 


18 


5.0 


7.19 


20.33 


37.35 


2.54 



Table A.2. Table 4. Same as Tab lA.ll except the WR2 mass loss 
rate is used. 



MwR [Mo] 


Mcc [Mo] 


Pdh] 


Pilh] 


Pilh] 


8 


1.4 


4.83 


13.67 




12 


1.4 


3.81 


10.77 




16 


1.4 


3.26 


9.21 




8 


3.0 


6.54 


20.02 




12 


3.0 


5.27 


15.76 




16 


3.0 


4.56 


13.48 




8 


5.0 


7.77 


25.84 


2.75 


12 


5.0 


6.39 


20.35 


2.26 


16 


5.0 


5.60 


17.41 


1.98 



either the strong mass loss of the WR star (WRa) or tidal spin- 
up (WRb + WRc). The system remains synchronized throughout 
its evolution till RLOF starts. 



Appendix A: Grid models setup 

Table lA.ll gives an overview of our grid of models using the 
WRl mass loss rate, the models using the WR2 rate are shown 
in Table |A3] 



Appendix B: Angular momentum evolution 

The detailed specific angular momentum profiles of the WR stars 
for each model (WRa, WRb and WRc) can be seen in Figures 
IB. II IB. 21 and IB. 31 These all show that during core helium burn- 
ing the specific angular momentum profile inside the star does 
not change significantly. Also clear is that the specific angular 
momentum of model WRa, does not increase during the evolu- 
tion of the system, while for models WRb and WRc the spin-up 
clearly works, i.e. the specific angular momentum increases. 

Fi gure s |BT4l and |B3] sho w the evolution of the orbital angular 
momentum and the degree of synchronization ( a»wR / Worb )■ 
The angular momentum of the orbit always decreases, due to 
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Fig.B.l. Specific angular momentum as a function of mass 
coordinate for our Cyg X-3 model WRa (1.4 Mq companion). 
Profiles are for the start of core helium burning (solid line), 50% 
He-depletion in the core (dotted line), 25% He-depletion in the 
core (dashed line) and the last model calculated when carbon is 
depleted in the core (dot-dashed Une). 



Fig.B.3. Specific angular momentum as a function of the 
mass coordinate for Cyg X-3 model WRc (1.4 Mq companion). 
Profiles are given for the start of core helium burning (solid line), 
50% He-depletion in the core (dotted line), 25% He-depletion in 
the core (dashed line) and the last model calculated when RLOF 
started (dot-dashed line). 




Fig. B.2. Specific angular momentum as a function of mass co- 
ordinate for Cyg X-3 model WRb (1.4 companion). Profiles 
are for the start of core helium burning (solid line), 50% He- 
depletion in the core (dotted line), 25% He-depletion in the 
core(dashed line) and the last model calculated when RLOF 
started (dot-dashed line). 



Fig. B.4. Orbital and stellar angular momentum as a function of 
time for our three Cyg X-3 models. The orbital angular momen- 
tum (solid lines) is plotted for model WRa (blue, solid), model 
WRb (green,dashed) and model WRc (red,dotted). The stellar 
angular momentum (dotted lines) has the same coloring for the 
different models. 
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Fig. B.5. The ratio of angular velocity and orbital angular ve- 
locity (jjwr / <j^bin for model WRa (solid blue line), model WRb 
(dotted green line) and model WRc (dashed red line). 



